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CO ' Abstract. In this paper, we prove second derivative estimates together with classical solvability for 

the Dirichlet problem of certain Monge- Ampere type equations under sharp hypotheses. In particular 
we assume that the matrix function in the augmented Hessian is regular in the sense used by Trudinger 
and Wang in their study of global regularity in optimal transportation [21] as well as the existence of 



o 
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J-H I a smooth subsolution. The latter hypothesis replaces a barrier condition also used in their work. The 

applications to optimal transportation and prescribed Jacobian equations are also indicated. 

1. Introduction 

<^ ; This paper is concerned with existence of globally smooth solutions to the Dirichlet problem for 
certain Monge- Ampere type equations. These equations arise in many applications, notably in optimal 
transportation, geometric optics and conformal geometry. The Dirichlet problem of the Monge- Ampere 
type equations under consideration has the following general form 

det{D^u- A{x,Du)} = B{x,Du), xeQ, 
u = ^p{x), X £ dQ, 



in 



(1.1) 



^ where Q is a bounded domain with smooth boundary in n dimensional Euclidean space M", is a smooth 
function on dQ. ^ is a given n x n symmetric matrix function defined on O x M"', and -B is a scalar 
i7~j valued function defined on x M". We use {x,p) to define points in x M" so that A{x,p) G M" x R"-, 
B{x,p) £ M. Also Du and D'^u denote the gradient vector and Hessian matrix respectively of a function 
ff^ u € C^(r2) . The equation in (|1.1|) reduces to the standard Monge-Ampere equation when ^ = 0. A 
solution u € C'^{Q) of (jl.ip is called an elliptic solution whenever D^u — A{x,Du) > 0, which implies 
^ B>0. 

^ . These equations have attracted significant interest in recent years; for a recent survey and the earlier 
• ^ ' history see, for example, [26]. The Heinz-Lewy example in [19] shows that there is no regularity 
^ ■ for the equation in (jl.ip without restrictions on A. When A = 0, the existence of globally smooth 
I solutions to the Dirichlet problem (jl.ip for smooth uniformly convex domains was finally settled by 
Caffarelli, Nirenberg and Spruck [l]|2j, Krylov [13] and Ivochkina \12\ I13j. For application to optimal 
transportation with quadratic cost function, the natural boundary condition is the prescription of the 
gradient image 0,* = Du(i}) which is equivalent to an oblique boundary condition for elliptic solutions. 
In this case local regularity for convex fi* was proved by Caffarelli [3] and the corresponding global 
regularity for smooth uniformly convex domains and target domains fi* proved by Delanoe [6j, 
Caffarelli ^ and Urbas [29j. The first breakthrough for general optimal transportation problems was 
obtained by Ma, Trudinger and Wang [18] who obtained an a priori interior second order estimate 
under an analytical structure condition on the matrix A, which they called A3, together with local 
regularity results under a generalized target convexity condition. These were subsequently extended 
to global estimates and regularity in [2H] under a weak form of this condition called A3w; (see PB] 
for further discussion). Following this Loeper [T7] found a geometric characterization enabling him to 
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adapt the proof of necessity of target convexity in [18] to prove that A3w is necessary for regularity in 
optimal transportation. 

In this paper we show that the condition A3w also suffices for global regularity in the Dirichlet 
problem not just in the optimal transportation case but also for general equations of the form (jl.ip . 
In this generality it is convenient to adopt the terminology from [23j and call the matrix function A 
regular if 

(1-2) Aj,kii^^P)^i^jVkVi ^ 0, 

for all {x,p) enx W\ ^,rj £ W\ ^±7], where Aij^ki = D^.p.^ij- 

We now formulate the main results of this paper. The first result is a global bound for second 
derivatives of equation p.l|) which replaces the barrier condition in |28j by the existence of a subsolution. 

Theorem 1.1. Let u € C^(Jl) n C^(fi) he an elliptic solution of problem in Q., where B G 

C^(il X M"), mi B > and A G C^(0 x M") is regular. Suppose also there exists a subsolution 
u G C'^{Q.) of equation Then we have the estimate 

(1.3) sup|L>^n| ^ C(l + sup[Z)2ni), 

n an 

where the constant C depends on n,A,B,n,u and sup {\u\ + \Du\). 

n 

From Theorem II. H we can infer a global bound for solutions of the Dirichlet problem (jl.ip . For this 
we adopt the approach from [71 18] and instead of explicit geometric assumptions on the boundary dO, 
as formulated for example in (2^, we assume more generally the existence of a subsolution with the 
given boundary trace. 

Theorem 1.2. In addition to the assumptions in Theorem \1. 11 suppose the subsolution u = ipon dQ 
with if G C^(i7),Or2 G C^. Then any elliptic solution u G C'^{Q) PI C^(0)) of the Dirichlet problem 
U.l\) satisfies the global a priori estimate 

(1.4) sup \D'^u\ < C, 

where the constant C depends on n, A, B,Q,,u and sup {\u\ + \Du\). 

n 

Existence theorems follow from Theorem ll.2l by standard methods. In general we need supplementary 
conditions on the matrix function A to control gradients and upper bounds of solutions. These can be 
reduced in the optimal transportation case, where the matrix A is generated by a cost function, namely 
a function c G C^(M"' x M") satisfying the conditions: 

(Al) For each x £ Q., p G M", there exists a unique y = Y{x,p) such that Dxc{x, y) = p; 
(A2) det Z?2 yC^ 0, for ah x G = Y{x,p),p G M". 

The matrix A is then given by 

(1.5) A{x,p) = Dlc{x,Y{x,p)). 

We formulate here a particular theorem for the general case which also embraces the main examples 
in optimal transportation. For this we assume a structure condition: 

(1.6) A{x,p)-^ -^lo{l + \p\^)I, 
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for all X G r2,p G and some positive constant /xq ^ind that the maximum eigenvalue of A{x,0) is 
non-negative. The last hypothesis can be replaced by a weaker condition that there exists a bounded 
viscosity supersolution to equation (jl.ip . 

We then have the following existence theorem. 

Theorem 1.3. Under the assumptions in Theorem together with the above structure conditions, 
there exists a unique elliptic classical solution u € C^{Cl) of the Dirichlet problem with u. 

The paper is organized as follows: In Section 2 we prove a fundamental lemma about barriers for the 
linearized operator and use it to establish the global second derivative bound in Theorem ll.il We also 
prove an analogous variant of the interior Pogorelev estimate of Liu and Trudinger in ^16j . In Section 
3, we derive the boundary estimates needed to complete the proof of Theorem 11.21 Our treatment 
of these boundary estimates also provides the proof of the corresponding results stated in [23] under 
the barrier condition introduced in p8] . In Section 4, we establish the necessary gradient and solution 
bounds for our existence theorems. Once we have the relevant a priori estimates, Theorem 1.3 follows 
immediately from the method of continuity in [11, Chapter 17]. Finally, we discuss separately the 
special cases of optimal transportation and prescribed Jacobian equations. 

2. POGORELOV TYPE ESTIMATES 

In this section we prove global and interior Pogorelov type estimates for second derivatives of equation 
(1.1). In particular we use the regularity of the matrix function A to construct a barrier function for the 
linearized operator from a subsolution. The resultant second derivative estimates then follow readily 
from the corresponding proofs in [28] and |16j . 

Let F[u] = log(det(D^u — A{x,Du))), it is known that -F is a concave operator with respect to 
D^u — A{x,Du). The linearized operator of F is defined by 

(2.1) L = F'^iDij - Dp^Ai.ix, Du)Dk), 

where F^^ = and {F^^} = {u;*-'} denotes the inverse matrix of {wij} = {uij — Aij{x, Du)}. Assume 

the function Aij{x,p) G C^ifl x M"), i,j = 1, • • • ,n, B{x,p) € C^{d x M"), B{x,p) = log{B{x,p)) and 
set 

(2.2) C = L-Bp^D,. 
We introduce a fundamental lemma here: 

Lemma 2.1. Under the above assumptions, suppose u is an elliptic solution of il.l]} . u is a strict 
elliptic subsolution of equation if A is regular, then 

(2.3) C (e^(^-")) ^ ei X] " 

i 

holds in 0, for positive constant K sufficiently large and uniform positive constants ei,C. 

Proof. Since n is a strict subsolution of (jl.ip . u satisfies 

F[u] = log(det(Z)^u - A{x, Du)) ^ log{B{x, Du) + Sq), 

for some 5o > 0. 

For any xq € il, let = u — ||x — xo|^. For e small enough, the perturbation function is still a 
strict subsolution and satisfies 

F[Ue] = log(det(L»2n, - A{x, Du^))) ^ log{B{x, Du^) + r), 
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for some positive constant r. 

Let V = u — u, Ve = — u. By calculation, we have 

Lv = L(ue) + - xoP) 

= eF** - eF^Wp^Aijix, Du){x - xo)fc + F'i{DijV, - Dp^Aij{x, Du)DkV,} 
^ ' ' = eF'' - eF^Wp^Aijix, Du){x - xo)k + F'^{Dij{u, - u) - [Aij{x, Du^) - Aij{x, Du)]] 

+F^i{Aij{x,Du;) - Aij{x,Du) - Dp^Aij{x, Du)DkVe}. 

By the concavity of F, we have 

F[u,] - F[u] ^ r^{Dij{u, -u)- [Aij{x,Du;) - Aij{x,Du)\}. 

By the Taylor expansion, for some 9 € (0, 1), we have 

Aij{x,Du^) - Aij{x,Du) - Dp^Aij{x,Du)DkVe 
= Dp^Aij{x,p)DkVe - Dp^Aij{x,Du)DkVe 

= ^Aij^kl{x,p)DkVeDlVe, 

where j5 = (1 - e)Du + ODy^, P = (1 - 9)Du + ODu^ and 9 G (0, 6). 
Thus, at X = xo, we have 

Lv ^ eF^^ + F[u,] - F[u\ + ^F'^ Aij^ki{x,p)DkvDiv 

= eF^^ + lF'iAj^ki{x,p)DkvDiv + log{B{x,Du) + T)-log{B{x,Du)) 
^ epii + ^F'^Aij^ki{x,p)DkvDiv - Ci, 

where Ci is a constant depends on B, Du, and Du. 

Let (j) = e^"" with positive constant K to be determined, we have 

= F'HDij^-Dp^Ai,{x,Du)Dk^) 
= Ke^'"Lv + K^e^'^F^WivDjV 

^ Ke^^{tF'' + ^^F'iAi.j^ki{x,p)DkvDiv - Ci + KF'WivDjv}. 

Without loss of generality, we assume that Dv = {Div,0, ■ ■ ■ ,0). Since A is regular, we have 

L(p ^ Ke^^ieF'' + ^F'^Aij^u{x,p){Div)'^ + KF^\Divf - Ci} 

^ Ke^''{eF'' + l ^ ^ F'^ Aij^ii{x,p){Divf + KF^\Divf - Ci}. 

i or j=l 

Since the matrix {F^^} is positive definite, any 2x2 diagonal minor has positive determinant. This 
implies 



Then the Cauchy inequality leads to the following inequality 

for positive constant rj. 
Thus, we have 

L(t> > Ke^^e^" - f ??F-|^ii,n(x,p)|(Dit;)2 - |-F"|.4i*,ii(^,p)I(^i^')' + KF^^D^vf - Ci]. 



obtain 



8»7 

Choosing ri small such that r? ^ rri -i — =vrr7T — vir and K large such that K ^ ^™ax 



L4) ^ Ke^^i^F" -Ci}. 



Thus, we have 



£ct) = L(l)- Bp^Di(t) ^ Ke^'"{-F^ - Ci} - Bp^D, 



which leads to the conclusion of Lemma 12.11 with ei = minQ{|i('e^^} and C = maxQ{CiKe^^ + 
Bp^D.cj)}. 



□ 



Remark 2.1. The function (f) = e^^- is global harrier function for the linearized operator C The 
inequality C (e^^-~"^) ^ ei ^ — C in Q will be used to control the second order global and interior 

i 

estimates in this section and the boundary estimates in the next section. 

Remark 2.2. It is a fairly standard calculation that a non-strict classical subsolution of a uniformly 
elliptic partial differential equation can be made strict using the linearized operator and the mean value 
theorem. For example, if u is a non-strict subsolution, then u-\- ae^^^ is a strict subsolution for small 
constant a and large constant b, (see [11], Chapter 3). This can also be done in a neighbourhood of the 
boundary, preserving boundary conditions by replacing xi by a defining function or more specifically the 
distance function. By using the mean value theorem we do not need concavity or convexity and the small 
constant a controls the uniform ellipticity near the subsolution. Hence we need only assume the existence 
of a non-strict subsolution in Lemma 2.1; the inequality 112. 3\) will still hold for the corresponding strict 
subsolution. Thus, the second order apriori estimates will also hold under the existence of a non-strict 
subsolution and we only need to assume a non-strict subsolution in the hypotheses for our theorems. 

Remark 2.3. In the previous paper [28], to obtain the second derivative bounds, a kind of global barrier 
condition, called A-boundedness in [23] is assumed, namely Q is called A-bounded with respect to u if 
there exists a function ip € C'^ipt) satisfying 

(2.5) [Aj^ - Dp^Ai,{-,Du)Dk^]Uj > lel', 

in 17, for all ^ € M". From Lemma 2.1 it follows that A-boundedness can be replaced by the existence 
of a subsolution for second derivative estimates. 

Remark 2.4. If F is substituted by another concave operator F[v\ = (det(-D^M — A{x,Du)))ri, the 
conclusion of this lemma still holds. The only difference is the form of F^^ = We have F^^ = 

^fw'^ if F[u] = {det{D'^u - A{x,Du)))^, while F'^ = w'^ if F[u] = log{D'^u - A{x,Du)). Here we 
make a convention that {w^^} denotes the inverse matrix of {wij} = {uij — Aij{x, Du)}. 

Now from Lemma 2.1, we obtain the global Pogorelov estimate, Theorem 1.1, by appropriate ad- 
justment to the proof of Theorem 3.1 in [55]. This is done by substituting the barrier function 
(f = exp(i^(u — u)) in place of the barrier from condition ()2.5p . That is we maximize the function, 

(2.6) e^p{^\Duf + b(t)}w^^, 

over and |.^| = 1, where w^^ = Wij^i^j = {uij — ^jj)CiCj, </> = e^^-~^^ with K as in Lemma 2.1, and a, 
b are positive constants to be determined. Then we arrive again at the estimates (3.11) and (3.13) in 
[28], but now with an additional dependence on n, and Theorem 1.1 follows. 

As remarked in the introduction we may similarly modify the proof of the interior Pogorelev estimate 
in [I6j to obtain the following variant. 
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Theorem 2.1. Let u € C'^{Vt) n C^'^{Vt) he an elliptic solution of problem M-l\) in Q., where B G 
(7^(0 X M"), iniB > and A € C^(17 x M") is regular. Suppose also there exists an elliptic suhsolution 
u € (7^(0) o/ equation and a degenerate elliptic supersolution uq € C"'^'^(r2) Pi C^''^{^) such that 

u = uq on dQ. Then we have for any il' (s O 

(2.7) sup|Z)2n|^C, 

n' 

where the constant C depends on n, A, B,Q,Q' ,u,UQ,and sup^duj + |I?n[). 

Note that if A{x,0) = 0, then our estimate (j2.7p agrees with the usual Pogorelov estimate [11], by 
taking uq = 0. 

3. Boundary estimates 

In order to complete the proof of Theorem II. 2| we need to obtain an a priori estimate for the 
Hessian D^u on the boundary dO, for solutions of the Dirichlet problem (jl.ip . First, we note from 
Remark 2.2 that we can assume the subsolution u is strict provided we restrict to a neighbourhood 
of dQ. As remarked in the previous section, this can be shown by adding a function ci exp(c2d(x)), 
where d{x) = dist{x,dil,) denotes the distance function in il, and ci and C2 are positive constants. 
Accordingly, we can retain Lemma 12.11 in a neighbourhood of dfl, which will suffice for boundary 
estimates. 

Next, we need to observe that the regularity of the matrix function A will be preserved under 
coordinate changes, that is, under a coordinate change, equation (jl.ip is transformed into an equation 
of the same form with transformed matrix A satisfying (jl.2p with respect to the transformed gradient 
variables. Specifically if we take a diffeomorphism, y = ip{x) with Jacobian matrix J = {ipij} = {Djipi}, 
then 

Wij{x) = ^kilPljDyi^y^U + Di{lljkj)Dy^U " Aij{x, JDyU). 

The transformed matrix function A' and scalar B' are thus given by 

\ A\j{-,p) = ^'Hi\Aki{'^-\Jp)-Dk{il^sl)Psi 
^ ' I B'{;p) = {det J)-^B{i;-\jp), 

where {^^•'} = J^^. The invariance of condition (1.2) follows readily from (3.1). 

Consequently, we may fix a point xq G dQ, which we take as the origin, and a neighbourhood J\f 
of Xq, such that T = M H dil lies in the hyperplane {xn = 0}, and the positive x„ axis points into $7. 
We then have by differentiation, Dai3u{x' ,0) = Dai3u{x' ,0), a, /? = 1, • • • , n — 1, x' = (xi, ■ ■ ■ ,x„_i), 
which leads to the double tangential derivative estimate \Da/3u{x)\ ^ C, a, /3 = 1, • • • , n — 1, on T. 

We then estimate the mixed tangential- normal derivatives DanU, ot = I,-- - ,n — 1, onTby dif- 
ferentiating the equation (jl.ip and using the barrier in Remark 2.1. With this barrier in hand, the 
derivation of these estimates is similar to the special case of the standard Monge- Ampere equation, as 
treated for example in [11]. For completeness and later reference, we carry out the details here. First 
by differentiation of the equation 

(3.2) F[u] = log(det(Z)\ - A{x, Du))) = B{x, Du), 
we have 

(3.3) F'-^iuaij - DaAij{x,Du) - Dpi^Aij{x,Du)uak) = Baix,Du) + Bpi^{x,Du)uo,k, a = !,■■■ ,n, 
which leads to 

(3.4) CDo,u = Baix,Du) + F'Wo,Aij{x,Du), 
where £ is defined by ()2.2p . Then we have for tangential derivatives, 
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(3.5) \C{Da{u-u))\^C{l + '^F''), a = 1,-- in 0. 

i 

We also observe that, on d^l close to the origin, the tangential derivative satisfies Da{u — u) = 0, 
a = 1, 2, • • • , n — 1, i.e., \Da{u — u)\ ^ C|x'p. This implies 

\Da{u — u)\ ^ C\x\'^ , on dil,s, 

where = O n ^^(O) is a smaU neighborhood of the origin, dCls = {dQ D Bs{0)) U{nn dBs{0)). We 
always choose 6 so small that d^l PI .65(0) C T. 

Let V = 1 — (j), where (j) = e^^-~^^ is the barrier function in Remark 2.1, then we have 

(3.6) Cv -ei^F'' + C, in n, and v = 0, on dn, 

i 

for positive constants ei and C. We consider a function of the form 

(3.7) = V + ^Xn - kXn^, 

where n and k are positive constants to be determined. We then have 

(3.8) ^ -j(l + ^F**), in ns, and ^^0, on 50^, 

i 

for k sufficiently large and fi, 6 sufficiently small. 

By modification of the function ip, we employ a new barrier function 

(3.9) i> = aip + b\x\^ , 

with positive a and b to be determined. By a direct computation, for a » 6 ^ 1, we have 

(3.10) £tp = aCi; + hC\x\^ ^ "(^ " + ^'')' 

i 

(3.11) |Z)„(u -u)| = ^ -0, onaJ]nB5(0), 
and 

(3.12) \Da{u-u)\<,h5'^ ^atp + hd"^ ='ip, on Vtr\ dBs{<d). 
Therefore, for a S> 6 S> 1 and 5-^1, there holds 

(^12,\ \CDa{u - u}\+ C-il) ^ 0, in Vts, 

\Da{u — u)\ ^ i/', on d^s. 

By the maximum principle, we have 

(3.14) \Da{u-u)\ ^ '4), in 1^^. 

At the origin, we have \Da{u — u)(0)| = ■0(0) = 0. For \Dci{u — u)|, taking x' = 0, dividing by x„ and 
letting x„ — )• 0, we get 

(3.15) |Z?«„(u-u)(0)| ^ aVn(O) ^ C. 

Thus, we obtain the mixed tangential-normal derivative estimate \Danu{x)\ ^ C, a = 1, ■ ■ ■ ,n — 1, on 
T. 

It therefore remains to estimate the double normal derivative DnnU and for this the regularity of the 
matrix function A is critical. Following the idea in |21j as used in [22], we fix a unit vector ^ G R"'"^ 
and set, (with respect to our transformed coordinates), 

W: = w[u\= 

(3.16) = [DapU- Aai3ix,Du)]^aCl3 

= [D^p^ - Aa,p{x, D'lp, Dnu)]iaii5 on T, 



with w > hy virtue of the ehipticity of u. Here D' = {Di, Dn-i) denotes the tangential gradient. 

An estimate from above for DnnU on T will follow from equation (jl.ip . provided we obtain a positive 
lower bound for w on T. To get this we note that for sufficiently large K the function w = w + K\x\'^ 
will take a minimum on T x S'^~^ at some point x € T, and for some unit vector ^. Extending (p by 
defining (p{x',Xn) = 93(2;', 0) for x„ > 0, it follows from the regularity of A, that the corresponding 
extension of w is a concave function of DnU. Consequently from the differentiated equation (3.4), for 
a = n, we have 

(3.17) ^ C7(1 + ^F"). 

i 

Now using the barrier in Remark 12. II as above, we obtain this time a one-sided estimate Dnw{x) ^ — C, 
that is 

(3.18) DnnU{x)Dp„Aai3{x, D'if, DnU)^aCl3 ^ C. 

To reach our desired estimate for DnnU, we then need to get a positive lower bound for the coefficient 
in (j3.18p . By the ellipticity of the subsolution u, we can fix a positive constant 8 for which 

(3.19) ^ 5, for aU x G T, |^| = 1. 
We also have 

"^["] ~ " -\^ap{x, D'lp, DnU) - Aap{x, D'lp, DnU)]CaCl3 

> -Dn{u - u){x)Dp,^Aal3{x,D'(f,DnU)^a^l3, 

since DnU > DnU from the strictness of u and again using the regularity of A. We also have an upper 
bound 

(3.21) D„iu-u)i^K, 
for a positive constant k. 

2' 



Combining (l3T9]) . ((3:20]) and (lOTD . we thus obtain for w[n] < ^ 



(3.22) i^p„Aa/3(S, i^V, i^n^^)eae"/3 > TT " 

Hence, we conclude from (j3.18p an estimate Dnnu{x) ^ C for a further constant C. Now utilizing 
equation (jl.ip again, we obtain an estimate from below for w{x) and finally an estimate from above for 
Dnnu{xo)- Since an estimate from below automatically follows from the ellipticity of u, we complete 
the estimation of D^u on 50 and the proof of Theorem 11.21 

Remark 3.1. In the above proof, we could have used, in place of i3.16\) . the function 

(3.23) w[u] = {detlwij]}"^ 

where indices i,j = l,--- ,n — 1, in direct accordance with the technique in [21] . 



Remark 3.2. In [23], the concept of domain A-convexity is introduced extending that of c- convexity in 
optimal transportation. In particular for A a given nxn symmetric matrix function defined on ^Ix M" 
we may define 17 to be uniformly A- convex, with respect to u if 

(3.24) [Daj + Dp^Aij{x, Du)-fk]TiTj ^ 60 

for all X £ d^l, unit outer normal 7 and unit tangent vector r and some positive constant 60. It then 
follows that ^} is uniformly A-convex if and only if, for any constant K > 0, there exists a defining 
function ip G C^(i7) satisfying ip = on dQ, together with the inequality Ii2.5\) in a neighbourhood of dU.. 



By using 93 in place of the barrier in Remark \2.1[ we conclude an estimate for the Hessian D'^u on dO, 
for solutions of the Dirichlet problem \1.1\ for uniformly A-convex domains, thereby proving Theorem 
2.2 in [23]. Note that the concept of A- convexity is also invariant with respect to coordinate changes. 
Note that the replacement of uniform A-convexity by the existence of a strict subsolution with the same 
boundary trace for second derivative estimates, as formulated in Theorem 2.1, is also pointed out in 

Remark 3.3. Theorem M.Si provides bounds for both D^u and the augmented matrix function w. By 
the equation M.l\j . we have the positive lower bounds for w in Ct. Thus, the uniform ellipticity of the 
operator F, with respect to u, follows easily from the positive upper and lower bounds of w. 

4. Existence theorem and some applications 

In this section, we complete the proof of the classical solvability result for the Dirichlet problem (11. ip 
and consider the applications to the optimal transportation and prescribed Jacobian equations. 

First we establish the necessary solution bounds and gradient bounds for Theorem 1.3. By the 
comparison principle, we have u ^ u in $7. The subsolution n is a lower bound of the solution. To obtain 
an upper bound, suppose u attains its maximum at xq G 0, that is Du{xo) = and D'^u{xq) ^ 0. By 
ellipticity, we then have A{xo,0) < 0, which contradicts the hypothesis that A{x,0) has a non-negative 
eigenvalue for each x € $7 and hence u must take its maximum on d^l. Thus we have the solution 
bound 

(4.1) \u\ ^ Kq, in Q, 

where Kq depends on ||m||loo(q) and ||</3||L°°(f7)- 

The latter condition is equivalent to constant functions being viscosity supersolutions, as defined in 
[20] . so more generally we could assume that there exists a viscosity supersolution u ^ </> on di^. 

We assume a quadratic bound from below for the matrix A to control the gradient of the solution, 
namely that A satisfy the following structure condition (jl.6p . 

A{x,p) ^ -/io(l + bP)/. 

Consider the function ip = e'^'^\Du\ for some k > 0. Suppose ^ attains its maximum at xq € Q, namely 
Di'ip{xo) = 0. Hence, at the point xq, we have DiuDiip = 0. By calculation, we have 

(4.2) k\Du\'^ + DijuDiuDjU = 0. 
From the ellipticity condition D^u — A[x, Du) > 0, we have 

(4.3) k\Du\^ + Aij{x, Du)DiuDjU < 0. 
By the structure condition (jl.6p of the matrix A, we have 

(4.4) kIDuI'^ - no{l + \Duf)\Du\^ < 0. 

Without loss of generality, we may assume Du{xq) ^ 1. The left hand side of (|4.4p will be nonnegative 
if we choose the positive constant n sufficiently large. This contradiction leads to the gradient estimate 

(4.5) sup \Du\ ^ C, 

where C depends on pL^, ||^i||Loo(n) and supg^ \Du\. 

On the boundary, the tangential derivatives of u are given by the boundary condition and the interior 
normal derivative bound from below is controlled by the subsolution u. The estimate from above follows 
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(4.7) 



from condition (jl.6p and the ellipticity, which impUes An > An, (see proof of Theorem 14.1 in [11]). 
Combining with (j4.5p . we have 

(4.6) \Du\ ^ Ki, in Cl, 

where Ki depends on |ju||ci(Q) a-^d ||</?||ci(n)- 

As we shall indicate below, for the special cases of prescribed Jacobian and optimal transportation 
equations, the above conditions can be relaxed. 

To complete the proof of Theorem 1.3, we have from Theorem 1.2 and the bounds (j4.ip and 
I, uniform estimates in (7^(0) for classical elliptic solutions of the Dirichlet problems: 

det{D'^u- A{x,Du)} = tB{x, Du) + {I - t) det{D'^u- A{x, Du)} in fl, 
u = if on dil., 

for ^ i ^ 1. Theorem 1.3 then follows from the Evans - Krylov second derivative Holder estimates of 
Evans, Krylov and Caffarelli-Nirenberg-Spruck, (see for example [llj . Theorems 17.26, 17.26'), and the 
method of continuity (|11). Theorem 17.8). The uniqueness assertion is immediate from the maximum 
principle. 

Remark 4.1. We remark that the classical solution u belongs to C°°(0) if both A and B are C°° 
functions on 0, x M", as well as d^l G C°°,(p £ C°°{dO,). Under the stated hypotheses it also follows 
from the linear theory [11] that the solution u G C'^'"(17) for all a < 1. In fact we have the stronger 
inclusion, u G T^^'^(O) for all p < oo and this would in fact suffice for the argument in Section 2. Note 
that we only use u £ C^{i),) in Section 3. 

Remark 4.2. As a special case for Monge-Ampere type equation (Ejp; existence results for the 
classical solutions to the Dirichlet problem for the equation det(uij + aij{x)) = ij){x) treated in the 
strictly convex domain in [U [2] can be generalized to non-convex domains. Also, the equation with 
general right hand side treated by Li [15j can be generalized to non-convex domains. 

Remark 4.3. For the relationship with equations arising in conformal geometry the reader is referred 
to [9] where a similar barrier argument for boundary estimates is used. 

In the remainder of this paper, we discuss the application to optimal transportation and prescribed 
Jacobian equations. Optimal transportation problems have received a lot of attention in recent years 
and for the basic theory, we refer to the the books, |301 l31j. Letting y be a mapping from x M x 
to M" and u a function in C^(0), we define a mapping T = Tu : O — t- M" by setting Tu = Y{-,u, Du). 
The corresponding prescribed Jacobian equation may be written the form 

(4.8) \dQlDT\ =■>!){■, u,Du) 

for a given positive function ^/^ on $7 x M x M". If det DpY ^ 0, we can then write equation ()4.8p . for 
elliptic solutions u, as an equation of Monge-Ampere type 

(4.9) de\.[D'^u - A{-,u, Du)] = B{-, u, Du) 

with matrix function A and scalar function B given by 

A{-,u,p) = -Yp-^{Y^ + Yu®p), 
^ ' B{-,u,p) = \detYp\-'^i^, 

so that when Y and ip are independent of u, (and B is positive), we obtain a Monge-Ampere type 
equation of the form (II. ip considered here. Optimal transportation equations are the special cases 
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where the mapping Y is generated by a cost function c satisfying conditions (Al) and (A2) as formulated 
in the introduction. We then have, in accordance with ()1.5p . 



(4.11) A{;p) = Dlci;Y), 

B{-,p) = \ det Dl yc\tlj > 0. 

The natural boundary condition for prescribed Jacobian equations is the prescription of the image 
r(f2). The existence of classical solutions is treated in [24j, which extends the optimal transportation 
case in [25]. However the classical Dirichlet problem in small balls has been used for local regularity 
arguments in [18 1 127 1 132]. Clearly, the estimates in Theorem ll.ll and ll.2l extend immediately to embrace 
prescribed Jacobian and optimal transportation equations under regularity of the matrix functions A 
given by (j4.10p and (j4.1ip . For the above solution and gradient estimates some relaxation of our 
additional hypotheses is possible. First we note that for any vector ^ G M", the quantity Tu ■ must 
assume its maximum and minimum values on the boundary dQ since det DT ^ 0. In the optimal 
transportation case this immediately provides a bound for the gradient in term of its boundary trace 
since Du = Cx{-,Tu) in Q. Also in the optimal transportation case, functions of the form u = c(-, yo) are 
solutions of the homogeneous equation so we automatically obtain solution bounds. For an arbitrary 
domain il, a condition of the type p.6p would still be required for the boundary gradient estimate or 
more generally the existence of a supersolution u satisfying u = ip on dQ. However if the solution is 
c-convex in 0,, that is at each xq £ ^l, there exists yo ^ 1^") such that 

(4.12) u{x) ^ it(xo) + c{x, yo) - c(xo, yo) 

in Q, then we have Tu{i^) C Tu{Q) by virtual of the monotonicity property (^18j, Lemma 4.3), and 
the global gradient bound follows immediately. From [27], we know that an elliptic solution will be 
c-convex in O if is c-convex with respect to each y E M"', that is the images Cy{-,y){Q) are convex in 
M". Consequently, we have the following existence theorem for optimal transportation equations. 

Theorem 4.1. Let c G C^(M" x M") be a cost function satisfying (Al), (A2) with regular matrix 
function A given by il.5\) and B G C^(i7 x M") satisfying mlB > 0. Suppose there exists a subsolution 
u G C'^{Cl) of equation in Q., satisfying u = Lpon dQ. Then if there also exists a supersolution 

u, = (f on dVl, or if i7 is c-convex with respect to all y G M", there exists a unique classical solution 
u G C^(0) of the Dirichlet problem with u ^ u in ^. 

Remark 4.4. A concept of generalized solution for optimal transportation equations, extending that 
of Aleksandrov for the standard Monge-Ampere equation, is also introduced in |18] and the Dirichlet 
problem for generalized solutions is considered in jlOj . Without the regularity assumption, this notion 
can be non-local as shown by Loeper [17j and thus is not really a weak form of the classical notion. 

Remark 4.5. These results also extend to the more general Monge-Ampere type equations |^.g| ) for A 
and B satisfying 

for all {x,u,p) G x M x M". In this generality. Theorem \4-l\ also extends to embrace prescribed 

Jacobian equations with mappings Y determined by generating functions as in ^25j . 
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